Nonsingular and accelerated expanding universe from effective Yang-Mills theory 
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The energy-momentum tensor coming from one-parameter effective Yang-Mills theory is here used 
to describe the matter-energy content of the homogeneous and isotropic Friedmann cosmology in its 
early stages. The behavior of all solutions is examined. Particularly, it is shown that only solutions 
corresponding to an open model do allow the universe to evolve into an accelerated expansion. 
This result appears as a possible mechanism for an inflationary phase produced by a vector held. 
Further, depending on the value of some parameters characterizing the system, the resulting models 
are classifled as singular or nonsingular. 
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I. INTRODUCTION 

The Friedmann-Lemaitre cosmological model described by the Robertson- Walker geometry with classical 

electrodynamics as its source, leads to a singularity at a finite time in the past It is usually identified as the origin 
I , of the universe. In fact, this singular behavior points out that, around the very beginning, the spacetime curvature 
becomes arbitrarily large, thus being beyond the domain of applicability of the model. It is expected that a quantum 
theory for gravity would circumvent the appearance of this curvature singularity by changing the predictions of the 
general relativity in the limit of large curvature. 

There are many proposals of cosmological solutions without a primordial singularity. Mostly they are based on 
, the fact that a collapsing phase can achieve a minimum, then evolving into a expanding phase. These nonsingular 

■ solutions are usually called bouncing cosmology. Such models are based on a variety of distinct mechanisms, such as 
\^ , cosmological constant [1| , non- minimal couplings Q , nonlinear Lagrangians involving quadratic terms in the curvature 
(N ■ i-ES, non-equilibrium thermodynamics ful] and loop quantum cosmology [T^, among others [IMl. In general, the 

I singularity theorems are circumvented by the appearance of a high (but nevertheless finite) negative pressure in 

■ the early phase of the universe. An initial singularity can also be avoided by assuming a quantum creation of a small 
Q\ I but finite universe, and hence a beginning of time [13, [2l|. For a recent review on bouncing cosmology see [2^ . 

, The issue of inflationary cosmology has also been largely considered in the literature. For a recent review on 

■ theory and observations, see 23| and references therein. Data from the Wilkinson Microwave Anisotropy Probe 
, (WMAP) observations [23i [25| have imposed several restrictions to the possible models describing an inflationary 

■ cosmology. Particularly, evidence [1^ for primordial non-Gaussianity in the temperature anisotropy of the cosmic 
' microwave background radiation (CMBR) seems to disfavor canonical single-field slow-roll inflation. Further, the 

spatial curvature of the universe is found to be small (with an spatial curvature parameter Hk of order 10^^), but 
not so small as predicted by the conventional inflationary models, where fix ^ 10^^. WMAP data seem to favor 
cosmological models with flat (Euclidean) or open (negative curvature) spatial sections. If the spatial curvature is non 
null, even being small, the implications for the time evolution of the universe can really be important, as described 
by general relativity. Recently [27| . inflationary cosmology and late-time accelerated expansion of the universe were 
considered in the context of non-minimally coupled effective Yang-Mills (Y-M) fields. It was shown that Y-M fields 
with a non-minimal gravitational coupling can produce an accelerated expanding phase in a flat (Euclidean spatial 
section) universe model. A complete analysis of the consequences for cosmology of the minimal coupling between 
gravity and effective Y-M fields has not been considered so far in the literature. The issue of inflation appears also 
in the context of nonsingular quantum cosmological models [1^ , where the inflationary phase appears as a quantum 
cosmological effect, instead of being produced by an inflaton field. 

In this paper, the cosmological consequences of the minimal coupling between gravity and effective Y-M theory 
are investigated. The energy-momentum tensor coming from the effective Y-M fields is assumed to be the dominant 
source of the matter-energy content of the early universe, which is assumed to be homogeneous and isotropic. The 
behavior of all possible solutions is examined for the three distinct topologies, corresponding to flat, closed (positive 
curvature) and open (negative curvature) spatial sections. It is shown that when Euclidean or closed sections are 
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considered, the solutions presenting expansion are decelerated and evolve into a static configuration. From this static 
phase the system can evolve into a collapsing phase if fluctuations on the scale factor are allowed. On the other hand, 
solutions corresponding to an open model do allow the universe to evolve into an accelerated expansion, thus providing 
a simple mechanism for an early inflationary phase produced by a vector field. Models presenting an inflationary phase 
produced by vector fields have been called as dark energy models ^29-31] . A graceful exit from the inflationary phase 
appears naturally because of a relationship between the scale factor and the magnitude of the Y-M field. Further, 
depending on the value of some parameters that characterize the system, the solutions turn out to be singular or 
nonsingular. Solutions exhibiting collapsing phases are also examined. 

The paper is organized as follows: the following section presents some general aspects of the homogeneous and 
isotropic cosmology, including the main equations. In Sec. IIIIl the energy-momentum tensor for general one-parameter 
Lagrangian density is derived. It is also presented the procedure of field averaging which renders the system a perfect 
fluid conflguration. Section IIVI deals with effective Y-M theory. In this context, the equation governing the time 
evolution of the scale factor is derived. In Sec. |Vl cosmological models resulting from the minimal coupling between 
gravity and effective Y-M theory are examined. There, a suitable functional form for the effective coupling is assumed. 
In Sec. IVIl we use the results obtained before in order to discuss the regime of small coupling with large mean flelds, 
for which the effective coupling presents a well known form. Finally, concluding remarks are presented in Sec. I VIII 

Heaviside non-rationalized units are used. Latin indices run in the range (1,2,3) and Greek indices run in the range 
(0, 1, 2, 3). Throughout the text the units oi c ~ 1 ^ h are used, unless otherwise stated. 



II. FRIEDMANN-LEMAITRE-ROBERTSON- WALKER UNIVERSE 



Maxwell electrodynamics as source of gravity in the homogeneous and isotropic Friedmann-Lemaitre cosmology 
l|, |3 leads to singular universe models. In this framework, this is a direct consequence of the singularity theorems 
19||, and follows from energy conservation law and Raychaudhuri equation [s^l • Let us assume the homogeneous and 
isotropic Friedmann-Lemaitre model as described by the Robertson- Walker metric |3| in the form 

ds^ = dt^ - (^^^ [dr' + ide' + sin^ 9 d^^)] , (1) 

where e = —1, 0, -1-1 hold for the open, flat (or Euclidean) and closed sections, respectively. Here {r,9,ip) are 
dimensionless comoving coordinates. In what follows the above mentioned framework will be simply referred as 
FLRW(60|. The 3-dimensional surface of homogeneity t = constant is orthogonal to a fundamental class of observers 
determined by a four- velocity vector field u'^ = Sq- 

For a perfect fluid with energy density p and pressure p, the energy conservation law and the Raychaudhuri equation 
read, respectively, 

p + 3ip + p)j=0, (2) 

^--^(p + 3p), (3) 

in which k is the Einstein gravitational constant {k = SirG) and 'dot' denotes Lie derivative respective to v'^, that is 
d/dt. It is worth to mention here that p + 3p> states the strong energy-dominance condition (SEC). For instance, 
in the case of A{t) describing an expanding phase, satisfying SEC means that A < 0. In this case we say that gravity 
decelerates the expansion. On the other hand, if this condition is violated the corresponding model will present a 
phase of accelerated expansion, as required for an inflationary cosmology. This issue will be further addressed when 
cosmological models are examined. 
Multiplying Eq. (O by AA we obtain 

AA + ^AA{p + 3p) = 0. (4) 



But AA{p + 3p) = —{A^p)\ where Eg. Q was used. Thus, introducing this result in Eq. and after integrating 
it, we obtain the Friedmann equation [6l| 
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which consists of a first integral of Eqs. ([2|) and If p is a known function of A{t), the above equation will determine 
the solutions for the scale factor. In other words, the time evolution of the scale factor is determined by Einstein's 
equations in terms of the matter-energy content of the universe. 



III. ENERGY-MOMENTUM TENSOR FOR ONE-PARAMETER NONLINEAR LAGRANGIAN 

This section provides a description of the tensor characterizing the matter-energy contents of the system, which is 
assumed to be described by a general one-parameter Lagrangian density. We begin by introducing some well known 
objects, which will be useful to set the notation. The strength tensor field F^^^ and the gauge field A^j^^ are related 

by 

F^^'J = ~ 9,4") + C^^^vlW , (6) 

where C"'^'^ represents the structure constant for a compact Lie group G. This tensor field can be conveniently defined 
in terms of the non-Abelian electric and magnetic sj^ ^ color fields as 

Fl.i> = v^.E\-) - v.E^;^ - rj^.'^^v^Bp . (7) 

In order to alleviate the notation the 'color' indices, above indicated by upper brackets, will be omitted in what 
follows. 

Let us assume the gauge invariant Lagrangian density to be a general function of the Lorentz invariant F 
F^'^F^^ = 2{B^ — E^) as L = L(F). The energy-momentum tensor for this class of one-parameter theories can be 
presented as 

T^i, = —ALpFfj a F°'„ — Lg^^, (8) 

in which Lp represents the derivative of the Lagrangian density with respect to the invariant F. 

Since the spatial sections of FLRW geometry are isotropic, we may consider that the color fields can generate such 
universe only if an averaging procedure is performed |34| . For this propose the system is assumed to satisfy the 
following requirements: (i) the volumetric spatial average of the color field strength does not depend on directions; 
(ii) it is equally probable that the products E''E^ , B^B^ and E^B^ (with i ^ j), at any time, take positive or negative 
values; (iii) there is no net flow of energy as measured by comoving observers. The above mentioned volumetric spatial 
average of an arbitrary quantity X for a given instant of time t is defined as 

(X)^ lim / X^/^dV, (9) 



V^Vo V 

with V ~ J y/—g cfix^, and Vq stands for the time dependent volume of the whole space. Similar average procedures 
have already been considered in [13, [33 - l36| . 

In terms of the electric and magnetic color fields, these requirements imply that 

{E,)= , (B,)=0, (10) 
{E^B,) = - ^-{E-B)g,,, (11) 

{E,E,) = - \E^g,„ (12) 

{B^B,) = - \B^g^,, (13) 

where we have defined E'^ = -E'E, and = -B'B,. Note that Eq. (HH) implies in {Ei Bj - Bi Ej) = 0. 

Applying the above average procedure to the energy-momentum tensor we obtain the following non null components: 

(Too) = -{ALpE^ + L)goo, (14) 
"4 ^ 

" 9^j, (15) 



Lf{E' - 2B') + L 



which can be presented as a perfect fiuid configuration with energy density p and pressure p as 

{Ty^v) ^ {p + p)Vf,Vu-pgy,i,, (16) 
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where we identify: 

p = -ALfE'^-L, (17) 
P = ^{E^ -2B^)Lf + L. (18) 

For the particular case of Maxwell electrodynamics, with L — —F/A, we obtain the classical result p = 3p — 
{E"^ + B'^)/2 (here E and B stand for the Abelian electric and magnetic fields). The fact that both the energy 
density and the pressure are positive definite for all times yields the singular nature of FLRW universes. The 
Einstein equations for the above energy-momentum configuration lead to ^38i] the classical solution for the scale 
factor A{t) — {A^t — ei^)^/^, where Ao is an arbitrary constant. 

For the case of nonlinear spin-one fields described by the Lagrangian density L = —F/A + aF^ + /3G^, which 
encompass the first order terms coming from one-loop QED, it can be shown [14| that, in the absence of electric 
field, p = (B2/2)(1 - 8aB'^) and p = {B'^/6){1 - AOaB^). This toy model presents solutions exhibiting a nonsingular 
behavior for the scale factor. Regards must be taken in applying such model to the description of the early universe. 
Higher order terms coming from the field expansion in the one-loop effective Lagrangian for QED may not be negligible 
in the regime of large mean fields, as it would be expected to occur in the early universe. Nevertheless, it could be 
considered in a specific phase where the above terms in the Lagrangian density would correspond to the dominant 
terms in the effective QED Lagrangian. 

Let us return to the general case determined by Eq. (jl6p . Taking in consideration the averaging procedure on 
fields, the energy conservation law stated by Eq. ^ reduces to 

^^E^§-(B^-E^) 
Lp at 

+ {B^ + -B')!- In \{B^ + E^)A^] ^ 0. (19) 
ot 

Summing up, the behavior of the scale factor A{t) in the FLRW cosmology dominated by a perfect fiuid given by 
Eq. (HH) is determined by Eqs. © and (fTQ]) . 

The dominant matter-energy contents of the universe in a very early phase is expected to be a plasma of quark 
and gluons (QGP). Such a phase would occur when the temperatures exceeded the value Aqcd ~ 200MeV [47|. At 
this regime the strong coupling is small and most quark and gluons only interact weakly. Soft modes could also exist 
but they would constitute only a small fraction of the total energy density. It is not clear by now what would be the 
better model to describe the evolution of the spacetime in this phase. The production of QGP is expected to occur in 
high energy collisions in the great accelerators, as in RHIC and LHC experiments (in BNL and CERN, respectively). 
The results to come from these experiments may shed some light in the understanding of the very early universe. In 
the following sections we shall examine a working model in which the Y-M effective Lagrangian density for quantum 
chromodynamics with one parameter background field will be considered to describe the matter-energy contents of 
the isotropic and homogeneous FLRW universe. The use of effective Y-M theory in the context of QGP and gluon 
plasma has recently been considered in [36|, . 



IV. ONE-PARAMETER EFFECTIVE LAGRANGIAN FOR YANG-MILLS FIELDS 

The effective Lagrangian density for quantum chromodynamics (QCD) in terms of the parameter background field 
F can be presented ,39. - .41.] in the form: 

\ F F 
Lym = --:^i-r^, 7 = log— (20) 

where the effective coupling 3(7) is implicitly given by means of 

rah) 1 

with /3{g) the Callan-Symanzik /3-function, p is the renormalization mass and g the gauge field coupling constant 
appearing in the basic QCD Lagrangian. This effective Lagrangian density should be taken as a classical model 
that incorporates several features of the quantum problem. It gives a sufficiently correct description of the quantum 
vacuum, opening a window to the examination of physically important configurations (40j . 
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In fact, there are many invariants of the Yang- Mills fields, depending on the specific gauge group [42[. The ansatz 
used to derive this effective Lagrangian density takes in consideration only the algebraic invariant F and imposes 
consistency with the trace anomaly for the energy- momentum tensor . 

Before analyzing the cosmological implications of this model, some remarks are in order. The study of the energy 
density for the effective action associated with this theory, revels that > would lead to a metastability of the 
vacuum. The interpretation for this behavior can be presented as follows: if a region in the system develops a large 
E field, it will quickly decay into a configuration where > E"^ |40j- Furthermore, we limit our considerations to 
the case in which only the average of the squared magnetic color field B survives [il, m Si-liil. This is formally 
equivalent to put E^ = in Eq. ([T2|) . With this assumption, the energy conservation law [cf. Eq. ([T9|) ] leads to 

Bit) = (22) 

Now, introducing this result in Eq. ([5]) we obtain the following equation governing the time evolution of A{t): 

The equilibrium solutions of this equation can be obtained by means of 



kB^ 1 

^1 



6 A^,^ - - ^''^ 
which applies for each possible value of the parameter e (0, ±1), corresponding to different topologies of the spacetime. 

V. COSMOLOGICAL MODELS 

In terms of an action principle formulation this work deals with a minimal coupling between gravity and Y-M fields, 
whose total action is given by 



5= / d 



j4. 



•:\f-det[g^]{LR + Lym), (25) 



where det[gp,^] is the determinant of the matrix whose elements are the components of the spacetime metric, given 
by Eq. ([!]), Lfi represents the well known Einstein- Hilbert Lagrangian density [IJ, and Lym is the effective Y-M 
Lagrangian density presented in Eq. ^U^. 

Now, let us examine the solutions for the time evolution of the scale factor A{t) for the case where the effective 
coupling g presents the form 

4 = 1 + (26) 
9 

As it will be shown in the next section, this form for the effective coupling correspond to the first two terms coming 
from the regime of small coupling with large mean fields [s^, 110, 13 ' where the constant bo is identified with a beta- 
function coefficient. The first term, the Maxwell-likc term, is really negligible at the mentioned regime, but it will be 
maintained here for future reference. Now Eq. (1231) reads 



a , b 



1/2 



with 



3 ' 



(28) 



b = ^ei/^o. (29) 

The derivative of Eq. (f27| . which provides an expression for the acceleration field associated to the function A{t), is 
given by (for A 0): 
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It is worth to mention that Eq. ([27)) can be presented as a second order differential equation in a more suitable form 
as j/ + 2a/y + e = 0, with y = A^{t). 

Formal solutions of Eq. (j27l) are implicitly given by 

/■^(*Va b \-^'^ 
t-U = ± In — - e dz, (31) 

Ja{u) \z z j 

which does not contain the equilibrium solutions. For each value of the parameter e, real solutions can exist only if 

^ln-^-6>0. (32) 

In the above condition, the equality gives the equilibrium solutions (or equilibrium points) of Eq. ([27| . In fact, at 
the equilibrium points we have 

e-^/" = ^. (33) 

As the scale factor A appears squared in this expression, all solutions will be of the form ztA^. Nevertheless, the 
negative solutions shall not be considered, since they do not correspond to physical results for the scale factor. 
Summing up, solutions for the time evolution of A{t) are given by Eqs. (pij) and (|33p . 

In what follows the graphic study of these solutions will be performed for each given spacetime topology. We are 
not assuming any particular normalization for the scale factor. Particularly, A = 1 does not imply present time. 



A. Euclidean section (e = 0) 



In the case of an Euclidean section, there is one positive equilibrium solution, which is given by Aq = It can 
be seen directly from Eq. (jSS]) by setting e = 0. Now, let us study the behavior of A{t) in the neighborhood of Aq. 
From Eq. ([27l) one can see that A{t) is real only \i A < 6^/^, otherwise it will be complex. Furthermore, in the region 
where A is real, it will be positive for the positive root of Eq. (I?7|) and negative for its negative root. In other words, 
the equilibrium solution Aq behaves as an attractor for the positive root (in the sense that any nearby solution tend 
forwards it) and as a source for the negative root (any nearby solution tend away from it). 

The slope field(6^] for the positive root of Eq. (P7)) is presented in Fig. [TJ This figure shows that any initial 




0.5 1 
t 



FIG. 1: Slope field for the positive root of Eq. (|27p in the case of e = (Euclidean section). The two solid curves correspond to the 
equilibrium solution Aq — fe^^^ (the right line) and a numerical solution of Eq. I|27p with an initial condition satisfying A{ti) < Aq. 
Any solution of the differential equation with initial condition satisfying A{ti) < b^^"^ evolves towards the equilibrium solution, which 
behaves as an attractor in the plot. There is no solution for which A > Aq. The expanding solution corresponds to a singular and 
decelerated expanding cosmological model. We set a = 1 and fe = 1. 



condition for A{t) satisfying A{ti) < Ao leads to a solution presenting a primordial singularity. Singularity is here 
understood in the sense that, looking backwards in time, A{t) achieves the value zero in a finite interval of time. 
After this singular origin, the scale factor increases in a decelerated rate towards the equilibrium solution Aq. The 
two curves plotted in this figure correspond to the equilibrium solution (the right line) and a numerical solution of 
Eq. (|27|) with an initial condition satisfying A{ti) < Aq. 
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FIG. 2: Slope field for the negative root of Eq. (|27|) in the case of e = (Euclidean section). The two solid curves correspond to the 
equilibrium solution Aq = (the right line) and a numerical solution of Eq. (|27p with an initial condition satisfying A{ti) < Aq. 
Any solution of the differential equation with initial condition satisfying A{ti) < Aq decreases quickly to zero. In this case the 
equilibrium solution behaves as a source. The contracting solution corresponds to a singular and accelerated collapsing cosmological 
model. We set a = 1 and b = 1. 

On the other hand, if the negative root of Eq. (P7|) is considered, the conclusions are quite different, as shown in 
Fig. [21 In this case the scale factor begins with a finite value, which can not be greater than Aq. Then, it decreases 
in an accelerated rate to the singularity. This solution corresponds to a collapsing singular universe model. In the 
figure, the equilibrium solution appears as the right line. The other curve corresponds to a numerical solution of Eq. 
(P7)) with an initial condition satisfying A{t.i) < Aq. 

As the equilibrium solution behaves as a source in the case of the negative root, it is unstable under small pertur- 
bation. Thus, if the system is initially in this state, any small fluctuation in the value of A(t) would lead the system 
to evolve into a collapsing phase. In this sense, the two roots could also be considered together. In this case, an 
initially singular universe would expand in a decelerated rate until the equilibrium solution Aq is achieved. Thus, it 
could evolve into a singular collapsing phase. 

As one can see, the equilibrium point Aq corresponds to an upper limit for the values the scale factor can take. 
The value of Aq is associated with the parameters appearing in the Lagrangian density and depends on the specific 
model for the effective coupling constant. There is no solution presenting an accelerated expansion in the Euclidean 
section model. 

Since Friedmann equation is invariant under time reversal[63|, solutions for the negative root of Eq. ((3T|) . as those 
appearing in Fig. [51 can be obtained from the solutions for the positive root simply taking f — — t in Fig. [1] In what 
follows only solutions corresponding to the positive root will be considered in plots. 

It was recently shown 27] that a non-minimal coupling of gravity with effective Y-M theory can produce accelerated 
expansion in the flat model. In fact, it was claimed earlier, in jSOl] . that even for the case of a minimal coupling, it 
would be possible to obtain an accelerated expansion by considering a primordial electric condensate. Nevertheless, 
following the reasonings of a condensate with only an electric color field (or even with > B^) would imply 
in vacuum metastability. Within the framework considered here, our results confirm that, for a magnetic condensate, 
there is no solution presenting an inflationary phase in the Euclidean model, if the minimal coupling is assumed. 

Finally, in the case of Euclidean section, the solution presented by Eq. (piT) can also be presented in terms of the 
inverse error function as 




(34) 



where Co is an integration constant and Erf (z) is the error function [5l| 




(35) 



Here we must add again the equilibrium solution A(t) — 6^/^. The results obtained for this case (Euclidean section) 
can also be obtained from the above prescription, i.e., by means of Eqs. (jM)) and (|35l) . 
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B. Closed section (e = +1) 

In this case the positive equihbrium solution of Eq. ()27|) is given by 

A+i = ^aW{h/a), (36) 

where W{z) represents the Lambert function 52]. Here the behavior of the solutions are quite similar to the case 
examined before, for the Euclidean section. The slope field for the positive root of Eq. (P7|) is presented in Fig. [31 
As occurs in the case of the Euclidean section, the equilibrium solution corresponds to an attractor in the case of 




0.1 0.2 0.3 0.4 

t 



FIG. 3: Slope field and two solutions for the positive root of Eq. I|27p in the case of e = 1 (closed section). The two solid curves 
correspond to the equilibrium solution A+\ (the right line) and a numerical solution of Eq. I|27p with an initial condition satisfying 
A{to) < A+i. Any solution of the differential equation with initial condition satisfying A{to) < A+i evolves into the equilibrium 
solution. There is no real solution for which A > A+i. The expanding solution corresponds to a singular and decelerated expanding 
cosmological model. We set a — 1 and 6=1. 

the positive root and to a source in the case of the negative root. In both cases no solutions can be found for which 
A > A^i. Figure |3] also presents the equilibrium solution (the right line) and a numerical solution corresponding to 
an expanding model. The latter is given by Eq. pip with a given initial condition satisfying A{ti) < A^i. This 
solution presents a singular origin, for which A{t()) — 0, evolving into a decelerated expansion in the direction of the 
static configuration determined by A+i. 

On the other hand, by considering the negative root of Eq. ((27)) . any given initial condition satisfying A{ti) < A+i 
leads to an accelerated collapsing model. 

The only difference between the solutions for Euclidean and closed sections consists on the value of the equilibrium 
point. Again, there is no solution corresponding to an accelerated expanding model. For the case of the positive 
root of Eq. (P?]). any expanding solution will evolve into a static configuration with A — A+i. Nevertheless, since 
the equilibrium solution is unstable for the negative root, the system in the static configuration could evolve into a 
collapsing phase, provided small fluctuations of the scale factor are allowed to occur. 

C. Open section (e = —1) 

In the case of an open section, the equilibrium solutions of Eq. (j27p are given by 

A_i = iy/aW{-b/a), (37) 

which is real only if (b/a) < (1/e). The mathematical constant "e" denotes the base of the natural logarithm. There 
are three distinct cases to be analyzed here, depending on the values of the parameters a and b. If a < eb no 
equilibrium solutions can be found; ii a = eb there will be only one positive equilibrium solution; and finally, if a > e6 
there will be two positive equilibrium solutions. We remember that these parameters are related with the physical 
quantities coming from the Lagrangian densities. These three branches can be clearly understood with the help of 
Fig. (|4]). In this figure the solid curve represents the function b/A'^ while the dashed curves represent the function 
exp(— A^/a) for different values of the parameter a. Equilibrium solutions are determined by the interception between 
the solid and the dashed curves [cf. Eq. ([33| ]. These curves are denoted by /, for a > eb; II, for a = eb; and ///, for 
a < eb. As one can see, if a > e6 there will be two positive and distinct values of A for which these curves intercept. 
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f (A) 




FIG. 4: The solid curve corresponds to the plot of the function f{A) — b/A^. The three dashed curves I, II and III, correspond 
to the plots of the function f{A) = e"'^ for a > eb, a = eb, and a < eb, respectively. The intercepting points between the solid 
and the dashed curves occur at the equilibrium solutions, which are given by Eq. p7[). 

In the limiting situation where the parameter a attains the exact value eb, there will be only one positive value of A 
for which the curves intercept to each other. On the other hand if a < e6 no coincident point can be found. 

Corresponding to each one of the above mentioned domains, the solutions of the differential equation (|27p . given 
by Eq. pip, will exhibit different behaviors. Therefore, it is worth to examine each case separately. 

1. a < eb: no equilibrium solutions 

li a < eb Eq. (j37p does not present any real solution and, as a consequence, there is no equilibrium point. Taking 
the positive root of Eq. (j27|) one can see, from the slope field in Fig. [SJ that any given initial condition implies in a 
singular solution. That is, for a given instant of time to the scale factor attains the singular value Alto) = 0. Looking 
backward, this means that the scale factor achieves the value zero in a finite interval of time. In this sense to can be 
thought as an initial time. For t > to the scale factor is an increasing function of t. From Eqs. p7| and (l30l) one 
can see that this expansion is initially decelerated, when A{t) < (e6)^/^, later evolving into an accelerated expansion 
when A{t) > (e&)^/^. The transition from the decelerated to the accelerated phase occurs at a certain transition time 
tT given by A{tT) = (e6)^/^, which depends on the parameters characterizing the effective Y-M theory. The solid 
curve in this figure corresponds to a numerical solution of Eq. (j27p and describes the above discussed expanding 
cosmological model with a primordial singularity. 

4 
3 

< 2 



1 
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FIG. 5: Slope field for the positive root of Eq. (|27p in the case of e = — 1 (open section) with no equilibrium solutions (a < eb). The 
solid curve represents a numerical solution of Eq. (|27p with a given initial condition and corresponds to an expanding cosmological 
model with two distinct phases: a decelerated expansion evolving into an accelerated expansion. We set a = e/2 and b — 1. 

If the negative root of Eq. (P7|) is take into account the opposite behavior is found. The possible solutions correspond 
to singular contracting cosmological models with a decelerated collapsing phase evolving into a singular accelerated 
collapsing phase. 
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2. a = e6; one equilibrium solution 

In this case Eq. (|37|) presents one positive equilibrium solution, which is here denoted by An. The slope field for 
the positive root of Eq. (P7|) is presented in Fig. [51 Depending on the initial condition, the system will behave quite 
differently. For an initial condition satisfying A{ti) < An, the corresponding solution will be singular [A(io) — for an 
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FIG. 6: Slope field for the positive root of Eq. (|27p in the case of e = —1 (open section) with only one equilibrium solution (a — eb). 
Three numerical solutions of this equation is also plotted. The right line represents the equilibrium solution An. Initial conditions 
such that A > An lead to solutions describing a nonsingular accelerated expanding universe. We set a — e and b = 1. 

initial time to], evolving into a decelerated expansion towards the equilibrium solution An (the right line in the figure). 
A numerical solution of Eq. ^17} with a given initial condition satisfying A < Ajj is presented in the bottom of Fig. [S] 
(the solid curve bellow the equilibrium solution) and describes a singular cosmological model presenting a decelerated 
expansion. On the other hand, if A{ti) > Ajj the corresponding solution will be nonsingular [A{t) > Ajj V t] and 
presenting an accelerated expansion. This cosmological model corresponds to a nonsingular accelerated expanding 
universe. For a given initial condition this solution is numerically obtained from Eq. ()27|) and is plotted in Fig. |6l 
appearing as the solid curve above the equilibrium solution (the right line) . These two kinds of solution are separated 
by the equilibrium solution, which behaves as an attractor for solutions satisfying A{ti) < An and as a source for 
solutions satisfying A{ti) > An. The equilibrium solution itself corresponds to a static cosmological model. As the 
equilibrium solution is unstable for those solutions satisfying the initial condition A > Ajj, a singular decelerated 
expanding phase could evolve into an accelerated expanding phase, provided small fiuctuations of the scale factor are 
allowed when the system is nearby the equilibrium point at A = An. 

Now, if the negative root of Eq. (P?]) is considered, the possible solutions will describe collapsing or static cosmo- 
logical models. As before, these solutions can be obtained directly by taking i — > — i in the above solutions. 

The solutions coming from the positive and negative roots can also be combined. For instance a collapsing phase 
can evolve into an expanding accelerated phase producing a bounce. Note however that this combinations do not 
correspond to a superposition of solutions. It can occur if small fluctuations on the scale factor are allowed to occur 
when the system achieves its equilibrium point. In this sense, the system initially in the equilibrium state can evolve 
into an accelerated expansion if a positive fluctuation of A(t) occurs, or yet, it could evolve into a collapsing phase if 
a negative fluctuation occurs. 

3. a > eb: two equilibrium solution 

Finally, when a > eb there will be two positive equilibrium solutions given by Eq. (P7)) . which are denoted by Aj 
and A''j, with A" < Aj. We notice that Eq. (I?/]) presents no real solutions in the interval Aj < A < A''j. Furthermore, 
as the parameter a takes the limiting value eb the both equilibrium points get the same value Ajj, which corresponds 
to the case of only one equilibrium solution, analyzed before. 

The slope field for the positive root of Eq. (P?]) is presented in Fig. [T] As one can see, the equilibrium point 
behaves as a source for those solutions satisfying the initial condition A{ti) > A\, while A'^ behaves as an attractor 
for those solutions satisfying the initial condition A{ti) < Aj. For initial conditions satisfying A{t.i) < A'j, the 
corresponding solution will be singular [A(to) = for a initial time to] evolving into a decelerated expansion towards 
the equilibrium point A". It corresponds to a singular decelerated expanding cosmological model. By the other hand, 
if A{ti) > A\ the corresponding solution will be nonsingular, with A(to) = A\, and accelerated expanding. It describes 
a nonsingular accelerated expanding universe. Since the equilibrium solution A\ is unstable, an initially static phase 
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FIG. 7: Slope field for the positive root of Eq. (|27p in the case of e = —1 (open section). The top corresponds to the region where 
A{t) > A*'j, while the bottom corresponds to the region where A{t) < Aj. There is not real solution between these two limiting 
values. The equilibrium solutions A'j and are represented by the right lines (with A'f < Aj). The two solid curves correspond to 
numerical solutions of the positive root of Eq. I|27p for different initial conditions. Particularly, any initial condition satisfying A > A'} 
leads to a nonsingular accelerated expanding model. We set a — 2e and b = 1. 



could evolve into a expanding one if small fluctuations of the scale factor are allowed. These two kinds of solution 
are separated by a region < A < A''j, for which no solutions can be found. In Fig. [7] the equilibrium solutions A'j 
and A''j are represented by the right lines. Two numerical solutions of Eq. (|27)) are also presented. They correspond 
to different initial conditions. The curve in the bottom was produced by an initial condition with A{ti) being smaller 
than Aj, while the curve in the top was produced by an initial condition with A{ti) being greater than A^j. They 
describe, respectively, a singular decelerated expanding model and a nonsingular accelerated expanding model. 

Again, if the negative root of Eq. (127]) is considered, the possible solutions will describe collapsing or static 
cosmological models, and they can be obtained by taking t ^ —t in the above solutions. 

Depending on the equilibrium state the system lies, it can evolve into an accelerated expansion [from the positive root 
of Eq. (|3ip ] or into an accelerated collapsing phase [from the negative root of Eq. (|3T]) ]. provided small fluctuations 
on the scale factor are allowed. In this sense, the solutions coming from the positive and negative roots can also 
be combined. For instance, a decelerated expanding phase can evolve into an accelerated collapsing phase, or yet, a 
collapsing phase could evolve into a expanding accelerated phase, producing a nonsingular model presenting a bounce. 



VI. SMALL COUPLING WITH LARGE MEAN-FIELDS 

For the case of small coupling, the beta function appearing in Eq. (PT|) can be expanded as j53l-l5 

m^-lbog' + hg' + ---, (38) 

where bo and bi are the usual ^-function coefficients defined in one- and two-loop orders [1^. Now, if we take the 
limit of large mean fields {F/ ji^ » 1) we obtain, from Eq. ([2T|) . 

=6o7-2^1og7 + --- (39) 

It has been argued (57i | that this expansion may also give the leading two terms in the effective action for weak mean 
fields {F/fi"^ << 1), because the magnitude of the effective coupling in Eq. ([5^ is small in both regions. 
Introducing Eq. ^ in Eq. (gHl), we obtain (40l.l5i| 

1 F 

Lym « --6oFlog — . (40) 
4 n 

The /?- function coefficient 60 is known as the asymptotic freedom constant, and can be presented as bo = 
(l/87r2)(ll/3)C2(G). In SU(3) QCD with Nf massless fermion flavors it reduces to 
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As one can see, the effective Lagrangian density describing the regime of small coupling with large mean fields, 
given by Eq. (j40p . corresponds to a particular case of the Lagrangian density considered in the framework discussed 
in the last section. In this regime the parameters a and 6, appearing in Eq. (j27p . are given by 

a - (42) 

b = — ^, (43) 

which are obtained by neglecting the first term in the right-hand-side of Eq. (|26l) . 

It is important to note that Eq. (j40p does not require further terms in the approximation as the field becomes 
larger. As the fields become larger, better Eq. (|40)) holds as the dominant contribution to the effective Lagrangian 
density. 

Before closing this section, some comments on the relationship between the physical parameters bo, fi and Bq 
and the behavior of the solutions for A{t) are in order. First, for the cases of Euclidean and closed sections, the 
behavior of the solutions are not affected by the values these parameters can take. On the other hand, for the 
case of an open section, the values of these parameters strongly determine the behavior of A{t). Particularly, they 
also make infiuence on the corresponding model being singular or nonsingular. Let us briefiy discuss this point in 
terms of the magnitude of the primordial magnetic color field for the case of expanding solutions in the open section 
model. If Bq < (18e^)^/^/ (k6oM^) the corresponding model will be singular with an earlier decelerated expansion that 
evolves into an accelerated expansion phase. If Bq = (18e^)^/^/(K6oM^) the corresponding model can be singular or 
nonsingular, depending on the initial condition A{ti) be smaller or greater than the equilibrium point (e 6)1/2 ^ ^^^^ 
respectively. In the case of A{ti) < An the system will be singular, evolving into a decelerated expansion until the 
equilibrium point is achieved. From this point it can evolve into an accelerated expanding phase or into a collapsing 
phase (also accelerated) depending on possible fluctuations of the scale factor A{t). If no fluctuation is allowed the 
system will last in the static configuration given by An. On the other hand, if A{ti) > An the corresponding 
model will be nonsingular and accelerated expanding. Finally, if Bq > {ISe^y^^ / (nboiM^) the expanding solution can 
be decelerated or accelerated, depending on the initial conditions for A{t) be smaller than Aj or greater than A^j, 
respectively. The singular solution is decelerated expanding and can evolve into an accelerated collapsing phase or 
into a static configuration. The nonsingular models are accelerated expanding. Such nonsingular models attains its 
minimum A^in at the greatest root of {-'{fiboB^/3)W[{18y^^ /(nbon'^B^)]}^/^ . 



VII. FINAL REMARKS 

WMAP observations [13, HH] have brought several implications for cosmology. For instance the parameter related 
with the spatial curvature flK of the universe is shown to be small, presenting a negative mean value. For some 
reference values [13], it is shown that the combination of the present WMAP data with the Hubble Space Telescope 
data implies in fix = -0.014±0.017; the combination of WMAP with SNLS data implies in fix = -0.011±0.012; the 
combination of WMAP with SNGold data irnplies in Hk = — 0.023 ± 0.014. These results seem to favor both flat and 
open models. Further, it was pointed out in [2^ that the WMAP data contain evidence against the null-hypothesis of 
a primordial Gaussianity. By analyzing the bispectrum of the WMAP up to the maximum multipole Imax = 750 they 
[26] found 27 < /tvi, < 147, with 95% of confidence level. This result disfavor canonical single- field slow- roll inflation 
predictions of /nl = 0, and suggests that alternative models for an early inflationary phase should be considered. 

In this work the cosmological implications of the minimal coupling between gravity and effective Yang-Mills theory 
was examined. An average procedure on matter fields was adopted in order to provide a possible way to describe an 
isotropic and homogeneous FLRW cosmology. The behavior of all solutions for the time evolution of the scale factor 
A{t) were graphically examined and the main results can be summarized as follows. Models presenting Euclidean (fiat) 
or closed (positive curvature) spatial section can describe a decelerated expanding universe, an accelerated collapsing 
universe, or yet a static universe. In all cases, but the static one, these models are singular. Only models presenting 
an open spatial section (negative curvature) can support an accelerated expanding universe. Further, depending on 
the initial conditions and on the values of the parameters a and b that characterizes the system, the corresponding 
model can be singular or nonsingular. The expansion will be accelerated if A > 0, which occurs provided A^ > eb. 
For the case of a < c6 the resulting model is singular and evolves into a decelerated expansion until A achieves the 
value (e6)i/2 in a finite time. At this point A = 0. Then, it evolves into an accelerated expanding phase. If a = eb, 
exactly, there is one equilibrium point given by A'^ = eb. Solutions with initial conditions A{ti) > (eb)^/^ result in 
a nonsingular and accelerated expanding model, while solutions with initial conditions A{ti) < (eb)^/^ result in a 
singular model presenting a decelerated expansion. Finally, if a > e6 there will be two equilibrium solutions. The first 
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one A'} smaller than {ebY^'^ and the second one greater than (e6)^/^. Between them, there is no physical solution 
for A{t). In this case, for any initial condition satisfying A{ti) < Aj, the corresponding solution leads to a singular 
model presenting a decelerated expansion. Here the expansion occurs until A{t) achieves the static configuration 
at Aj. Since this lower equilibrium point is unstable for solutions coming from the negative root of Eq. (|27p . the 
system can evolve into an accelerated collapsing phase. Otherwise, if A{ti) > A^j the corresponding solution leads to 
a nonsingular cosmology presenting an accelerated expansion. In this case A{t) attains its minimum value at A\ in a 
time to- Solutions presenting collapsing phases can be obtained from the negative root of Eq. pi|) . and if combined 
with the expanding solution, provide a mechanism for a bounce. It is worth to stress that, considering the assumed 
framework, solutions presenting accelerated expansion only appear in an open universe model. Thus, it provides 
a possible mechanism for an early infiationary phase produced by a vector field. The solutions can be singular or 
nonsingular, depending on the initial conditions and also on the values that the physical parameters can take. 

Let us pick up a specific application of the results obtained in this manuscript. Let us take, for instance, a nonsingular 
expanding solution from the open-section model with the matter fields described by the Y-M Lagrangian density in 
the regime of small coupling with large mean fields [as described in Section IV C 3[ with the parameters a and b given 
by Eqs. ([^^ and (021) ]• In this case the universe would come from a previous decelerated contracting phase, which 
is a solution coming from the negative root of Eq. (I27p . and after achieving a minimum volume, determined by Aj, 
would expand in an accelerated rate. During this expanding phase the strong energy dominance condition is violated 
and, mathematically, the accelerated phase would last forever. Nevertheless, Eq. ((22|) shows that as A{t) increases, 
B{t) decreases. Thus, the system may achieve a phase in which the regime of large mean-fields does not apply. This 
fact shows that a mechanism for a graceful exit is supported by this framework. In fact, when F approaches the 
value /Lt^, the dominant term in the effective Lagrangian would be Maxwell- like. This is formally equivalent to neglect 
the second term in Eq. ([26l) . which would lead to a phase of decelerated expansion. In this sense, a inflationary 
phase of a nonsingular homogeneous and isotropic universe (63 | with a graceful exit appears naturally in the minimal 
coupled Yang-Mills fields in an open FLRW cosmology. Matematically, as consequence of considering Y-M matter 
field as the dominant energy contents, all the accelerated phases described in the open section models last forever. 
However, as suggested above, these solutions should be considered only as describing an early phase of the universe. 
As the universe expands the plasma is expected to evolve into another form of matter, which would dominate a next 
phase of the evolution. In this case, a complete scenario should consider the transition between this primordial phase 
to the next one, as radiation for example. In the case of a singular and accelerated expanding model (as presented 
in Section IV C 1[) the amount of infiation is a quantity of interest to be derived. This quantity can be obtained as 
N — \og[A{tend) / A{tbegin)] whcrc tend and tbegin represent the time coordinate at the end and at the begin of the 
inflation phase, respectively. Denoting the critical value of the magnetic field by Bcr, bellow which radiation would 
be the dominant form of matter, we obtain TV = {l/2)log{fj.'^ /V^eBcr)- 

Closing, few remarks are in order. First, the results obtained in this manuscript presenting infiationary solutions 
only occur in models with open spatial sectio n [65 1 ■ Infiationary models obtained from non- minimal Y-M theory have 
been considered only for flat spatial section [231 ■ One important feature of the solutions presenting an inflationary 
phase in the minimal coupling is that such phases can also appear in the context of nonsingular models. Thus, in such 
cases the amount of inflation is not a fundamental aspect to be considered in the solution of the initial conditions 
problems stated by standard cosmology. The motivation to work with models with non-flat spatial sections may be 
set by further observations from WMAP and PLANCK projects. It should be stressed that the spatial curvature has 
currently been found to be small but not necessarily zero, and the consequences of a non-zero curvature (even really 
small) can be relevant in the evolution of the universe, as shown here. Second, all the solutions presenting accelerated 
phases obtained in the context of miminal Y-M theory may be naturally obtained as limiting cases from non-minimal 
couplings with an open spatial section. Therefore, if non-flat models is considered in non-minimal coupling cases 
our results should be recovered in the corresponding regime. Finally, the anisotropics in CMB found in current 
data indicates that the early phase of the universe was not perfectly smooth. This fact must be connected with the 
inhomogeneities in the distribution of galaxies on large scales, as measured by Two Degree Field Galaxy Redshift 
Survey |59|. The very relationship between large scale structure and CMB is still an exciting issue in cosmology. In 
order to apply the results obtained here in the understanding of structure formation, the perturbations around the 
smooth background should be considered for each particular solution. In this context the power spectrum in the 
Fourier space would be the most important statistic. In the context of the solutions presented in this manuscript, 
two different situations could occur. If singular models are considered, inflation would be the responsible for the 
generation of scalar and tensor perturbations. By the other hand, if nonsingular models are considered, the primordial 
perturbations could possibly be produced in a previous contracting phase, thus evolving into the expanding phase by 
means of the bounce. The derivation of the perturbed Einstein equations in the context of minimal Y-M theory and 
its application in the study of the above mentioned aspects deserve future investigation. 
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